disp('fmdhlq');

n = 5;
t = 0.2;

y = fn_solve_basic(t, n);
disp(y);
disp(y - n*log(1+exp(-y)) - t);



function y = fn_solve_basic(RHS, n)

y0 = RHS;

while 1
    
   y = RHS + n*log(1 + exp(-y0));
   fprintf('y = %f \n', y);
   
   if abs(y - y0)<1e-12
       break;
   else
       y0 = y;
   end
    
end


end